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A preconditioning of the Cauchy-Rlemann equations which results In a 
second-order system is described. This system is shown to have a unique 
solution if the boundary conditions are chosen carefully. This choice of 
boundary condition enables the solution of the first-order system to be 
retrieved. A numerical solution of the preconditioned equations is obtained 
by the multigrid method. 
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Introduction 


In this paper we are concerned with the Cauchy-Rieraann equations. The 
two-dimensional form of these equations may he written as follows: 

U + AU = B in n, (la) 

X y 

U*n = g on r, (lb) 



// p = / g ds (2) 

n r 

which follows from Gauss' divergence theorem. This condition guarantees that 
the problem has a unique solution. We shall restrict ourselves to the case 
where ^ is the unit square. 

The Cauchy-Rlemann equations occur in many problems in fluid mechanics. 
The velocity components u and v are usually known as the primitive 
variables. It has been found that the numerical solution of the Cauchy- 
Rlemann equations can be enhanced by introducing a stream function, which 
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satisfies a Poisson equation. There are many fast, direct methods for solving 
such equations, see for example, Buzbee, Golub, and Nielson [2] and Hockney 
[6] . The emergence of multigrid methods means that there now exist fast and 
efficient iterative methods which are competitive with these direct methods, 
see for example, Foerster and Wltsch t^]. However, the primitive variable 
formulation does possess advantages since it allows treatment of more general 
flows. Ghll and Balgovlnd [5] have proposed a direct method to solve for the 
primitive variables in a problem where these variables are periodic in the 
x-directlon. 

This work was motivated by an examination of a least squares formulation 
of the Cauchy-Rlemann equations (Fix and Rose [8]). In two dimensions Gatskl, 
Grosch and Rose [9] employ an Iteration scheme due to Kaczmarz to solve the 
resulting problem but the method was found to be too inefficient to serve as a 
basis for a three-dimensional treatment of the incompressible Navler-Stokes 
equations. The method described here could provide the basis for a fast 
Cauchy -Riemann solver in three dimensions. 

This technique provides insight into the treatment of more difficult 
problems. Johnson [10] has successfully used a closely related method to 
solve the full Euler equations in straight channel with a 10% half-thick 
circular arc airfoil mounted on its lower wall. 

Here we describe a formulation of the problem which results in two 
Poisson equations, one for each of the primitive variables. Section 2 
describes the preconditioning and shows that with suitable extra boundary 
conditions the new problem possesses a unique solution which satisfies the 
Cauchy-Rlemann equations. Section 3 contains details of the solution 
procedure and describes the components of the multigrid method used. Section 
4 contains numerical results and computational details. 
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2. Preconditioned Equations 


We may write the system of equations (1) in the form 



We consider a preconditioning of this system obtained by preraultlplylng Eq. 

T 

(4) by the matrix D . This yields the following equation: 


(D*’^D)U = D^B 


(5) 


where 


and 



2 2 2 2 2 
The symbol V represents the Laplaclan operator 9 /9x + 9 /9y . This 

preconditioning has led to a system of second-order equations and so for the 

problem to be well-posed we require an additional boundary condition. In the 

remainder of this section we show that It is natural to take v - u =5 as 

X y 

the extra boundary condition. For simplicity we assume homogeneous boundary 
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data, i.e. U*n = 0 where n Is the unit outward drawn normal to the 
domain. There Is no loss of generality In doing so since U*n = 0 can always 
be achieved by an appropriate translation. In the following analysis we make 
the assumption that u, v, p and K have sufficiently many continuous 
derivatives . 

Theorem 1 

The solution of the system of equations (5) with the conditions that 
U*n = 0 and v^ - u = ? on the boundary of 5^ Is unique . 


Proof ; Let u^^^, v^^^ and u^^\ 

( 1 ) ( 2 ) , ( 1 ) - ( 2 ) , , 

u' ' ^ u and v * v . Let 

then <t> and satisfy the following 


v(2) 

<}> = u 


be two solutions of (5) with 


V 


2 


<!> = 0, 


'(' = 0 , 


In 0, 




X 



on r, 


= 0, on X = 0 and x = 1, 


'I' = 0, on y = 0 and y “ 1. 


The condition that = 0 along y = 0 and 

along these lines since \(; - <|) =0 there. 

X y 

problem for <j>: 


y = 1 Implies that = 0 
Thus we obtain the following 
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<)) = 0 



(|) = 0 


Figure 1 


It follows by a maximum principle (Protter and Weinberger [7]) that 
i}> = 0 In n. In a similar way it can be shown that 4* = 0 In 
the solution of the system (5) Is unique. 

Theorem 2 

The solution of the second-order system (5) satisfies Eg. (1) . 
Proof ; The following identifies hold: 

y2 _3£4.1i=3P_i2 

^ 3x 9y 3x 3y ’ 

^ 3y 3x 3y 3x * 

where P = u + v - P and Q = v - u - 5. 

X y X y 

and Q = 0 In . 


Thus 


We aim to show that P H 0 
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If U is a solution of the system (5) then 

r 

/ 

V^Q = 0 in On F we have the 

Q = 0 in by the uniqueness of the 
solution of Laplace's equation with Dlrichlet boundary conditions which 
follows from a maximum principle (see Protter and Weinberger [7]). 

Now 9P/8x = 0 and 9P/9y = 0 in implies that P is a constant, 
with value M say, in 

Using Green's theorem we obtain 


9x 


9y 


9y ^ 9x 


It follows, by differentiation, that 
boundary condition Q = 0. Therefore 


// (u + V - p)dx dy = / U*n ds + // P dxdy = 0. 
X y P 


( 6 ) 






The first term on the right-hand side of Eq. (6) is zero since U*n = 0 

around F while the second term is zero by the compatibility condition. 
Therefore, since u + v - P = M, a constant, we must have M = 0. 

» X y ’ 

Therefore, P = 0 in We have thus shown that the solution of the system 

(5) satisfies the Cauchy-Rleraann equations. 

It is certainly true that the solution of the Cauchy-Riemann equations 
satisfies the system (5) since the latter can be viewed as being obtained by 
differentiation of the former. The results of this section carry over to 


three dimensions and these will be reported elsewhere 
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3. Ifethod o£ Solution 

The system of equations (5) represents two Poisson equations for u 
and V. The problem has now been cast into a form, which is amenable to many 
of the standard Iterative methods of solution. We have chosen to use the 
multigrid method since such methods are fast and efficient. 

The domain 1s covered with a square grid of mesh length h = lA^ 

where N is a positive Integer. Each of the Poisson equations: 



l£ _ ii 

9x 3y * 


7^v = 


3p 


H 

9x 


are discretized using standard second-order central difference approximations. 
Along X “ 0 and x = 1 we have Dlrlchlet boundary conditions for u» Along 

the boundaries y = 0 and y = 1 the finite difference formula is 

constructed by eliminating exterior Imaginary points using a central 
difference approximation to the Neumann boundary condition ~ 

Since v is given along y = 0 and y = 1 we know the value of v^^ along 

these lines. Thus, we obtain a problem for u which does not contain any 

unknown values of v. Similarly, we obtain a problem for v which does not 
contain any unknown values of u. Therefore, the system for u and v 
becomes uncoupled and we can solve for the one independent of the other. 

We consider a multigrid method of solution to this problem using the 
Correction Storage algorithm of Brandt [1]. Let G, be a sequence of 
grids approximating the domain with corresponding mesh sizes hj^,***,h^. 
Let h|^ = ^^k+1 ^ “ 1,***,M-1. The problem is discretized on each 

grid Gj^ as described above. We describe below the components chosen in the 
multigrid procedure. 
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(1 ) Relaxation 

Pointwlse Gauss-Seldel relaxation Is used with the points ordered In the 
checkerboard (even-odd) manner. This Is one of the components of the "non- 
standard" multigrid techniques Introduced by Foerster, Stuben, and Trottenberg 
[3]. At the end of the relaxation sweep the residuals on any grid G|^ are 
zero at the odd points since at this stage all the odd-polnt equations are 
simultaneously satisfied. 

(2 ) Flne-to-coarse Transfer 

Checkerboard Gauss-Seldel produces highly oscillating residuals and so It 
Is Inadvisable to simply transfer the residuals by Injection to a coarser 
grid. Instead we transfer the residuals by half-weighting to the coarse grid. 
If the points of the coarse grid are also even points on the finer grid then 
the transfer of the residuals reduces to half- Inject Ion. 

(3) Coarse-to-flne Transfer 

Bilinear Interpolation Is used to transfer the correction to the fine 
grid to provide a new approximation there. 

We use a version of multigrid which begins with an approximation on the 
finest level of discretization and then cycles between grids until the 
convergence criterion has been attained. The decisions to switch grids are 
made automatically by Internal checks In the algorithm. The decision to 
switch to a coarser grid Is made If 

11^ ^'*■^11 2 > n 11-^^ II 2 , 

where n < 1 Is a positive constant and denotes the residuals of the 

discrete equations on grid Gj^. The decision to return to a finer grid Is 
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made according to the following criterion: 

Ik < 611^^11^ , 

2 2 ’ 

where 5 < 1 is a positive constant. 

4. Humerlcal Besults 

The raultigrid method defined in the previous section was used to solve 
three test problems. These three problems are characterized as follows: 

(i) P = 0, 5 = 0, 

u = X, V = -y; 

(il) p = 0, K = (m^ + n^) sln(nTTx) sln(mTry) , 

u = sln(nTTx) cos(mTTy), v = - (-^) cos(nTTx) sln(mTry); 

(iil) P = 0, ? = [ - n^TT^)/a]e sin(mrx) , 

u = e sln(mrx), v = (mr/a)e cos(nTrx). 

The first problem is homogeneous with a linear solution. The driving 
function, boundary conditions and solution of the second problem are 
oscillatory and this behavior Increases as m and n increase. The third 
problem presents a test of the multigrid method on a boundary layer problem. 

In all of the numerical experiments we took the finest grid to be such 
that N = 32 and considered a total of five grids in the raultigrid context. 
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The initial approximation on was taken to be zero everywhere except where 
the solution Is specified by the boundary conditions. With this initial 
approximation, the initial residual U^Qli 2 was computed and the convergence 
tolerance on was chosen to be 


e 


M 


10"^ IL^qII2 • 


We define a work unit to be the computational work in one relaxation sweep 
over the finest grid Gj^, The value of the switching parameters n and <S 
were chosen to be 0.6 and 0.2 respectively. 

The results of the algorithm applied to u for problems (1), (il), and 
(ill) are shown in Tables I, II, and III respectively. We give the number of 
work units required to attain the convergence criterion, the asymptotic 
convergence factor X, and the error in the numerical solution 

for u on Gj^. Similar results were obtained when the algorithm was applied 
to the primitive variable v in the above cases and indeed the number of work 
units was the same in all cases. We have used the notation a*b - c for 
a*b X 10“*^. 

It can be seen that the method exhibits the usual multigrid behavior by 
examining the as)rmptotic convergence factors obtained. The accuracy of the 
discrete approximation to problem (ii) decreases as m and n increase as 
one might expect since the number of mesh points per wavelength of the 
solution decreases. 

The approach described in this paper to precondition the Cauchy-Riemann 
equations together with the given choice of boundary conditions, which is 
based on firm theoretical results, has led to a problem which facilitates a 
fast numerical solution. This fast solution procedure has been obtained using 


a multigrid method 
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Table I. Details of method for problem (1) 





Work Units 

A 

llerrorll2 

13 

0.49 

0.45-4 


Table 

H 

H 

• 

Details of method 

for problem (11) 

m 

n 

Work Units 

A 

II error II 2 

1 

1 

14 

0.52 

0.43-3 

2 

2 

14 

0.51 

0.17-2 

3 

3 

14 

0.51 

0.37-2 

4 

4 

14 

0.51 

0.66-2 

5 

5 

12 

0.47 

0.10-1 

7 

7 

11 

0.43 

0.20-1 

3 

7 

12 

0.46 

0.14-1 


Table 

III. Details of method 

for problem 

(ill) 

m 

a 

Work Units 

T 

B error H 2 

1 

0.4 

11 

0.45 

0. 12-2 

1 

4.0 

11 

0.45 

0. 74-3 

1 

40.0 

12 

0.46 

0.92-2 
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